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Abstract. In the measurement of atmospheric v e and fluxes, the calculations 
of the Super Kamiokande group for the distinction between muon-like and electron- 
like events observed in the water Cerenkov detector have initially assumed a mis- 
identification probability of less than 1 % and later 2 % for the sub-GeV range. In the 
multi-GcV range, they compared only the observed behaviors of ring patterns of muon 
and electron events, and claimed a 3 % mis-identification. However, the expressions and 
the calculation method do not include the fluctuation properties due to the stochastic 
nature of the processes which determine the expected number of photoelectrons (p.e.) 
produced by muons and electrons. Our full Monte Carlo (MC) simulations including 
the fluctuations of photoelectron production show that the total mis-identification rate 
for electrons and muons should be larger than or equal to 20 % for sub-GeV region. 
Even in the multi-GeV region we expect a mis-identification rate of several % based on 
our MC simulations taking into account the ring patterns. The mis-identified events 
are mostly of muonic origin. 

Comparing the intensities per kiloton-year between the 33 kiloton-year (the first 
half of the total exposure of 79.3 kty) and 46.3 kton-year (the latter half) exposures, 
moreover, the respective values for e-like event numbers, /i-like event numbers and the 
ratio for (e-like number) /(/x-like number) are not consistent within the given statistical 
uncertainties in the sub-GeV and multi-GeV regions for the two subsamples. It should 
be also remarked that results for the whole zenith angle distributions of e-like and of 
/i-like events in the sub-GeV ranges and those of e-like and /i-like(FC) events in the 
multi-GeV ranges are inconsistent for the 33 kton-year and 46.3 kton-year exposures. 
We also studied the significance probability of the results after subdividing the whole 
zenith angle distribution into upward (— 1 < cos 6* < 0) and downward (0 < cos 9 < 1) 
going events. Their significance probabilities also support our claim that the method 
used by the Super K. group is inadequate for the identification of muons and electrons. 
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1. Introduction 

The most recent Super Kamiokande (Super K.) results for the study of atmospheric 
neutrinos and the neutrino anomaly before the large accident of the detector have been 
presented for the 79.3 kiloton-year (kty) exposure by Kajita and Totsuka |L] and at 
the 27th International Cosmic Ray Conference, Hamburg, August 2001 f2|. Using their 
reports, we have compared the respective intensities/kty as shown in Table 1 for the 
sub-GeV results (visible energy, E vis <1330 MeV) and for the multi-GeV results (E vis 
>1330 MeV). These results are subdivided into exposures of 25.5 kty(analysis A) and 
25.8 kty(analysis B) g g|, also 33.0 kty g, |, and 46.3 (= 79.3 - 33.0) kty 0, §, as 
given in the published papers and reports. 

Table [I] shows a comparison of various event samples for e-like event numbers, /i- 
like event numbers, (e-like+/i-like) numbers and the ratio for (e-like number) /(//-like 
number) in the sub-GeV and multi-GeV regions ||. The actual event numbers are 
shown in the second and the fourth column. Disagreements between the related two 
exposures are indicated by respective identifier marks. Earlier data for the sub-GeV 
region collected during the period from May 1996 to October 1997 have been analysed 
independently by group A for the exposure of 25.5 kty and group B for the exposure 
of 25.8 kty, respectively ||. The comparisons are shown in the upper part in table 
Both the event rates per kty for /i-like number, (e-like number + /i-like number) and the 
ratio for (e-like number) / (/z-like number) appear to be in conflict above the respective 
statistical uncertainties. Therefore, our main purpose for the present considerations is 
to compare those event rates for both exposures. The systematic uncertainties depend 
on the different kty values, but they are smaller than the statistical ones. 

In both A and B groups, the event rates for /i-like event numbers marked by < and 
for (e-like + /i-like) event numbers marked by > in the sub-GeV range do not agree. 
Also, the ratios for (e-like) / (/i-like) events marked by A are inconsistent in spite of the 
fact that the same data sample is analysed. However, the Super K. group finally has 
considered only the result of analysis A disregarding analysis B which has given larger 
number of /i-like events than analysis A. 

Various observed intensity values of the 33.0 kty exposure (including the data in 
25.5 kty run ||) are compared to those of the 46.3 kty in the central part of Table 1. 
The former is the first half in the total 79.3 kty exposure and the other is the latter half 
of the 79.3 kty exposure. The e-like intensities marked with O are n °t quite consistent 
with each other within the respective statistical uncertainties. 

Similar comparisons are given for the multi-GeV region. The /i-like "fully contained 
(FC) events" relate to all of the visible energy being contained in the inner detector and 
"partially contained (PC)" refers to events with some of the visible energy in the outer 
detector. For the identification procedure in the multi-GeV region, the Super K. group 
utilizes the characteristic Cerenkov ring pattern of muons and electrons. A /i-like event 
produces a clear outer ring edge of Cerenkov light while an electron-like event creates 
a fuzzy ring pattern due to the production of cascade showers and multiple scattering 
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Table 1. Comparisons of various event samples for sub-GeV and multi-GeV in the 
Super Kamiokande detector. 
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of low-energy electrons The comparison between the 33.0 kty and 46.3 kty cases 
are given in the lower part of Table 1. The intensities of e-like events marked with ® 
appear to exhibit a large difference (of the order of %) comparable to the sub-GeV case. 
Accordingly both intensities for (e-like events + /i-like events) and also the ratio (e-like 
intesity)/(/i-like intesity), marked with and © respectively, are different beyond their 
statistical uncertainties. However, intensities for /z-like (FC) and (PC) events for both 
exposures are consistent within the respective uncertainties. 

2. Study of the e/fi identification capability of the Super K. detector 

The only experimental test of the e//i identification capability was carried out by 
part of the Super K. group using the 12-GeV proton synchrotron at KEK, Japan p. 
For the test, they used a 1000 ton water Cerenkov detector. The water tank was 
cylindrical (about 10 m in height and 10 m in diameter) and its inner detector had 
a cylindrical volume of 9.6 m in diameter and 9.3 m in height. It was equipped with 
380 photomultiplier tubes (PMTs with 50-cm dia.). On the other hand the Super 
Kamiokande detector has a cylindrical volume, 39 m in diameter and 42 m in height, 
containing about 50 000 metric tons of water. The fiducial volume of 32 000 metric tons 
is cylindrical of 33.8 m in diameter and 36.2m in height. It is equipped with 11 146 
50-cm PMTs (1773 PMTs each on the upper and lower surfaces and 7600 PMTs on the 
sides) which cover 40 % of the inner surface. It is therefore conceivable that the large 
scale differences between the two detectors may lead a larger e//i mis-identifications for 
the full Super K. detector than the one that was obtained for the 1000 ton exposures. In 
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addition, the experimental data from the test set-up are of low statistics and limited only 
to sub-GeV events. There are no results for multi-GeV events because of the smallness 
of the 1000 ton detector. 

The Super K. group has performed a Monte Carlo (MC) calculation to establish a 
relation, Nmc{&, Pe), between the average number of photoelectrons (p.e.), the electron 
momentum P e , and the angle of incidence 9. The angle was between the particle 
direction and the center of circular area of 50 cm diameter. In the framework of this 
simulation it was assumed that the electron with a given momentum was started from 
a vertex at the center of sphere of 16.9 m radius. The relation N M q(9, P e ) is based only 
on the average number of p.e., neglecting the fluctuations around its average. By using 
this relation the expected p.e. number produced by the electron in the i-ih PMT on 
the inner surface of the cylindrical volume was calculated with the following expression 

(16 9 in\ ' ^ / l \ 

—j — j x exp f -j] x 

This expression was given as equation (6.7) in Kasuga's PhD paper [[K]]. In this 
expression, the exponent of 1.5 for (16.9 m/Zj) should have been 2.0. Also a p.e. number 
distribution for muons is given in a similar expression (6.8), which is 



N i:exp (direct) = (ct M x 1 — -g^— x sin 2 ^ + N iMock {Bi) ] 

x exp (--^j x /(e). 

The explanation for the denominator in the second case is incorrect in Figure 6.6 of 
Kasuga |10[ and Figure 2.14 of Sakai [11], afterwards being corrected by Ishihara in 
Figure 5.7 |12|| . Also they define the probability function suitable for the Particle 
Identification (P.I.), given by 



Prob(N exp , N obs ) = exp 



{N obs - N t 



\2 S 



exp/ 



2ttct 1 V 2a 2 )' 

with a 2 = 1.2 2 x N exp + (0.1 xN exp ) 2 . The first term with the factor 1.2 in a comes 
from the actual PMT charge resolution and the second factor takes into account the 
gain uncertainty in the PMT. The likelihood functions L e for e-like events and 
for yU-like events shown in their expression (6.12) are obtained with the probability 
Prob(N exp , N obs ) for studying the e//i identification capability of the events, L e = 
log(n 0i<( i. 5x0c ) Probi(e)) and L M = log(n 0i<( i. 5x0e ) Prob^n)), with an opening angle 
of 1.5 x6 c from the reconstructed particle direction. By adding information of the 
Cerenkov opening angles combined with the information of the ring pattern, the used 
method is improved with a PID (Particle Identification) parameter as given in their 
expression (6.16) to obtain a reduced mis-identification probability. Using the Monte 
Carlo simulation by applying these expressions to the experimental data obtained at 
the 12-GeV proton synchrotron, the group concluded |L0l |13| that they were able to 
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identify electrons and muons with mis-identification probabilities of less than a few % 



in the momentum range 250-1000 MeV/ c. In the similar study in Sakai's PhD paper [IT 



however, the mis-identifications for fi — > e are given as several % and for e — > // as a few 
% as shown in Table 7.1 and Table 7.2 in his paper using the same expression with the 
exponent 2.0, without considering stochastic processes of cascade showers. The Super 
K. group does not apply the ej \i identification method for multi-GeV region. They 
argue that muons in multi-GeV range can produce very large p.e. numbers compared 
to electrons of the same momentum, even increasing for larger momenta. Therefore, the 
e//i identification procedure applied for the sub-GeV range is unnecessary for multi-GeV 
events according to their estimates. 



3. Our full MC simulation 



The Super Kamiokande collaboration has determined the relations relevant for the 
mis-identification probabilities ( 6j,P e j or 8j,P^j) by only using the average number 
of photoelectrons. Of course, it is natural to expect that there are many events 
where the observed number of p.e. fluctuates around the average number. These 
fluctuations are not only caused by Poisson distributions, but also by the stochastic 
processes in the electromagnetic interactions of electrons and muons. Moreover, it is 
also an oversimplification to use only ionization and knock-on electron production for 
the expected p.e. numbers in the muon expression (equation (6.8) in the paper) ||10|| . 
However, in order to describe correctly the behaviors of fluctuations, additional energy- 
loss processes, like direct electron pair production and bremsstrahlung cannot be 
neglected. In the cylindrical fiducial volume, also, different path lengths of Cerenkov 
light due to different vertex positions and incident angles contribute to the fluctuations of 
the p.e. numbers. Therefore we simulated 100 events each with different vertex positions 
and determined the average p.e. number, while the p.e. number obtained by the Kasuga 
expression (6.8) depends only on the distance between U between the i-th PM and the 
vertex. To clearly demonstrate the effect we show the total photoelectron numbers 
for muons and electrons in the sub-GeV and multi-GeV ranges in their dependence on 
various vertex positions (different z-values) at x = y = cm and for vertical incidence 
in figure [I] and [| 

In view of such a complicated situation, we have performed a full MC simulation 
including all kinds of fluctuations using the EGS code for the Super Kamiokande 
detector. The EGS simulation can be initiated for electrons or photons which originate 
from a muon. As far as all other necessary quantities (water transparency, vertex 
position resolution and angular resolution etc.) are concerned, we follow the procedure 
used by the Super K. group. The step size of track lengths in our Monte Carlo is taken as 
0.028 radiation lengths for electron events. Also we use the probability function defined 
by Kasuga [[TIJ and determine the likelihood functions. As an appropriate coordinate 
system, we consider the x- and y-axes to be in a horizontal plane and the z-axis to 
extend perpendicularly downward. So the upper circular surface is at z= cm and the 
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electron momentum (GcV/c) 



Figure 1. Dependence of the average p.e. numbers on electron momenta in the 
sub-GeV and multi-GeV ranges at various z - values for vertical incidences. 



1: z= 200cm 



20J 




1 2 3 4 5 



muon momentum (GeV/c) 

Figure 2. Dependence of the average p.e. numbers on muon momenta in the sub- 
and multi-GeV ranges at various z - values for vertical incidences. 



lower one at z= 3400 cm. The non-linear behaviour in both figures reflect effects of 
vertex positions (corresponding to different z values), by which some tracks of events 
escape from the fiducial volume. 

As a reference, figures |3| and |] show the fluctuation pattern of electrons and of 
muons of 300 MeV/c for the vertex position of z= 200 cm and 6 = 0° obtained with 
1000 MC events. The average p.e. number for electrons with 300 MeV/c is larger than 
that of muons with the same momentum, but the situation is reversed at 5 GeV/c. 
The full black points in the respective figures show the Poisson distribution calculated 
for each average p.e. number where the track started at the respective vertex position. 



Identification problems of jj and e events in the Super K. 



7 
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Figure 3. Dependence of the p.e. number distribution for electrons of 300 MeV/c 
with a vertex position at z=200cm. 



average number of 
• p.e .= 569.3 




p.e. number 



Figure 4. Dependence of the p.e. number distribution for muons of 300 MeV/c with 
a vertex position at z=200cm. 



Figures [5] and ^ give the results for the same momentum at a different vertex point 
(z = 1700 cm). Figures and |] are for 5 GeV/c at z = 200 cm and figures || and 
|j"0| for 5 GeV/c at z = 1700 cm. In these figures we have shown only events for the 
minimum momentum value in sub-GeV range and the maximum momentum value in 
multi-GeV range. These results show that electrons of 300 MeV/c started at z = 200 
cm and 1700 cm show comparatively larger deviations from the Poisson distributions 
than those of muons under the same condition. On the other hand muons of 5 GeV/c 
at both given z-values exhibit larger fluctuations than electrons. Also, for the inclined 
incidence directions similar fluctuations occur with even stronger deviations compared 
to the vertical case. Accordingly, the fully simulated MC results clearly prove that 
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Figure 5. Dependence of the p.e. number distribution for electrons of 300 MeV/c 
with a vertex position at z= 1700cm. 

• • average number of 
event number »p e = 866 7 
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p.e. number 

Figure 6. Dependence of the p.e. number distribution for muons of 300 MeV/c with 
a vertex position at z=1700cm. 



the fluctuations owing to the stochastic processes can produce a serious impact on the 
identification problems of muon and electron events not only in the sub-GeV region but 
also in the multi-GeV region through their effect on the ring patterns. 

3.1. OurL^ — Le distributions 

We inspected the — L e distributions for electron and muon events with various 
momenta in the sub-GeV range and also some in the multi-GeV range for various starting 
positions (z- values), keeping x = y = 0. Then, the momenta of electrons and muons were 
selected such as to produce almost the same p.e. numbers. The mis-identification rates 
between both event types are obtained from their estimated results as obtained from the 
L M — L e distributions. For this purpose 500 muon and electron events each were used. 
As shown in table the corresponding mis-identification probabilities are summarized 
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Figure 7. Dependence of the p.e. number distribution for electrons of 5 GeV/c with 
a vertex position at z=200cm. 
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Figure 8. Dependence of the p.e. number distribution for muons of 5 GeV/c with a 
vertex position at z=200cm. 



and identified by various symbols. The mis-identification rates for muons (% values) as 
explained below the table are defined in the following way: [(muon numbers misidentified 
as electrons) — (electron numbers misidentified as muons)]/ (original number of muons). 
The mis-identification rates appear to be larger than the ones given in Kasuga's paper. 
From the results shown in the table, one might consider it difficult to reliably distinguish 
electron from muon events in particular in the sub GeV region, although the situation is 
improving with increasing momentum. However, these mis-identification rates are only 
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Figure 9. Dependence of the p.e. number distribution for electrons of 5 GeV/c with 
a vertex position at z= 1700cm. 
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Figure 10. Dependence of the p.e. number distribution for muons of 5 GeV/c with a 
vertex position at z= 1700cm. 
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Table 2. L^-L e distributions between electrons and muons in sub-GeV and some of 
multi-GeV. 
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the resultant mis- identification rate for muons; 

O ~ 5 %, o: 6 ~ 20 %, A: 21 ~ 60 %, <: 61 ~ 95 %, x: 96 ~ 100 % 
*: means L^-L e distribution with the opposite signs. 



As an example, we show in figure |TT] the distributions of — L e for 300 MeV/c 
electrons and for 500 MeV/c muons (the events with the chosen momenta will radiate 
almost similar p.e. numbers on average) for z = 2700 cm. In some cases it is difficult 



to distinguish betweem the two event types. Figure [12] gives the same quantity for 2000 
MeV/c electrons and for 2000 MeV/c muons at z = 1200 cm which allows to distinguish 
between the two event types. But, this is only true for vertical incidence keeping x 
= y = 0. One should notice that the — L e distribution for muons is much wider 
than that for electrons. This aspect is very likely caused by stronger fluctuations due 
to stochastic interaction processes of muons rather than fluctuations due to the shower 
development of electrons. As seen in table |2|, the — L e distributions of electrons and 
muons would lead to a mis-identification of < few % for muons in the region of z = 200 
cm ~ 2200 cm typically in the multi-GeV region. However, these results are only true 
for events started at the given vertex position of x = y = and for vertical incidence. 
For z >2200 cm, mis-identification rates of > several % for muons are obtained, and 
it has to be kept in mind that some muons could escape through the bottom surface 
of the detector. Moreover, such muons with oblique incidences can generate — L e 
distributions differently shaped compared to those for vertical incidence. 
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Figure 11. L M — L e for 300 MeV/c electrons and for 500 MeV/c muons at a starting 
point of x = y = and z = 2,700 cm. 

Number of events 




-100 100 200 300 400 

Figure 12. L^ — Le for 2000 MeV/c electrons and for 2000 MeV/c muons at a starting 
point of x — y = and z = 1,200 cm. 

3.2. Our results with PID parameter 

Following Kasuga's method of Particle Identification (PID) using the pattern and 
opening angle of the Cerenkov ring, furthermore, we calculated the PID parameter 
values. For the sub-GeV region, we obtained the % values for muon deficits in table 
[3] using a number of 9240 for electrons and 14 910 for muons starting at vertex points 
from z = 200 cm to 3000 cm with their incidence direction 6 = 0° along the z-axes (x = 
cm, and y = cm, 50 cm, 100 cm, 250 cm, 500 cm and 750 cm). The event numbers 
are taken from the intensities of the v e and energy spectra |15| . Then, according to 
our simulation, 4312 electrons and 6958 muons are estimated to stop up to z = 1400 
cm. e — > fi means the electron is misidentified as a muon, and /i — > e the opposite. The 
table shows substantial shifts from muons to electrons. 

For instance, in case of both events starting at z = 200 ~ 3000 cm and at x = y = 
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Table 3. Summarized PID results for distributions between electron and muon events 
in the sub-GeV region. 
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3485 


688 


4735 


increase 


+78.9 % 


-48.9 % 


+33.0 % 


-20.6 % 


+78.4 % 


-48.6 % 


+43.8 % 


-27.1 % 




7) x=0, y=1000 cm 


8) x=y=0 cm, 6 = 10° 


event number 


4312 


6958 


9240 


14910 


4312 


6958 


9240 


14910 


e — * /i, /x — * e 


25 


3307 


596 


4570 


6 


3995 


211 


5694 


increase 


+76.1 % 


-47.2 % 


+43.0 % 


-26.7 % 


+92.5 % 


-57.3 % 


+59.3 % 


-36.8 % 




9) x =y=0 cm, = 20° 


10) x=y=0 cm, 6 = 30° 


event number 


4312 


6958 


9240 


14910 


4312 


6958 


9240 


14910 


e — * /i, /x — * e 


3 


4197 


173 


6421 


1 


4226 


142 


6877 


increase 


+97.3 % 


-60.3 % 


+67.6 % 


-41.9 % 


+98.0 % 


-60.7 % 


+72.9 % 


-45.2 % 



with 9 = 0°, the percentage of misidentified muons is —7.1 % = (143 — 1202)/14 910 and 
+ 11.5 % = (—143 + 1202)/9240 for electrons. From these values, we obtain (number 
of muons misidentified as electrons) — (number of electrons misidentified as muons) a 
deficit of 7.1 % of muons and an 11.5 % increase of electrons, i.e. the interactions of 
muon neutrinos would be systematically misinterpreted as those of electron neutrinos. 
For reference, some results for inclined incidence (8 = 10°, 20° and 30°) at vertex 
positions of x = y = cm are added. Integrating numerically these values for muon 
deficits considering the number of neutrinos as given by the v e (z/ M ) energy spectrum, we 
obtain a total mis- identification between electron and muon events for starting positions 
at z = 200 ~ 3000 cm to be at least > 20 % for the sub-GeV region. 

In the multi-GeV region, we used 1897 Monte Carlo events for electrons and 2380 
MC events for muons as predicted by the respective neutrino fluxes. Similar PID 
parameters were calculated with vertex positions of z = 200 cm to 1400 cm and of 
z = 200 cm to 2900 cm. Some of the muons started at z > 1400 cm would escape 
from the fiducial volume. Table [| again indicates the deficits of muons similar as in the 
sub-GeV region. This leads us, thus, to conclude that even in the multi-GeV region the 
analysis procedure used by the Super Kamiokande group shifts about more than several 
% of muons into the electron category. 
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Table 4. Summarized PID results for distributions between electron and muon events 
in the multi-GeV region. 



1V1U.1L1 - Vjrc V 

region 


z=200 ~ 

1400 cm 


z=200 ~ 

2900 cm 


z=200 ~ 
1400 cm 


z=200 ~ 
2900 cm 


e 




e 




e 


M 


e 


M 




1) x=0, y=0 cm 


2) x=0, y=100 cm 


event number 


949 


1190 


1897 


2378 


949 


1190 


1897 


2378 


e —>■ /i, /i — > e 


5 


18 


95 


167 


23 


37 


88 


188 


increase 


+1.4 % 


-1.1 % 


+3.8 % 


-3.0 % 


+1.5 % 


-1.2 % 


+5.3 % 


-4.2 % 




3) x=0, y=250 cm 


4) x=0, y=500 cm 


event number 


949 


1190 


1897 


2378 


949 


1190 


1897 


2378 


e — > /x, — > e 


4 


41 


81 


73 


14 


57 


88 


236 


increase 


+3.9 % 


-3.1 % 


-0.4 % 


+3.0 % 


+4.5 % 


-3.6 % 


+7.8 % 


-6.2 % 




5) x=0, y=750 cm 


6) x=y=0 cm, 9 = 10° 


event number 


949 


1190 


1897 


2378 


949 


1190 


1897 


2378 


e — > jit, jit — •> e 


18 


118 


102 


309 


7 


13 


84 


58 


increase 


+ 10.5 % 


-8.4 % 


+ 10.9 % 


-8.7 % 


+0.6 % 


-0.5 % 


-1.4 % 


+1.1 % 




7) x=y=0 cm, 9 = 20° 


8) x=y=0 cm, 9 = 30° 


event number 


949 


1190 


1897 


2378 


949 


1190 


1897 


2378 


e [i, [i e 


31 


14 


148 


38 


5 


41 


81 


73 


increase 


+ 1.8 % 


-1.4 % 


-5.8 % 


+4.6 % 


+3.8 % 


-3.0 % 


-0.4 % 


+0.3 % 



For checking with the ring patterns, we inspected 300 events each randomly selected 
from 1897 electron events and 2380 muon events starting at vertex positions from z = 
200 cm to 1400 cm for vertical incidence along the z-axis (x = cm, and y = cm, 
100 cm, 250 cm, 500 cm and 750 cm). This resulted in a number of cases where it was 
difficult to distinguish muon from electron events which corresponded to 5 % and 11 %, 
respectively, in each 300 event sample. 



4. Comparisons of the zenith angle distributions in the sub-GeV and 
multi-GeV ranges 

Subdivided into the 33 kty (the first half) and 46.3 kty (the last half) expressions from 
the total 79.3 kty exposure, we compare both the observed zenith angle distributions for 
the sub-GeV and multi-GeV ranges using x 2 -tests [0. All the numbers of events can be 



read off the Super K. papers. Although the Super K. group uses 10 bins for the observed 
zenith angle distribution of the complete 79.3 kty exposure [[HJ, the 10 bins are only 
available for the 79.3 kty exposure for the sub-GeV and multi-GeV ranges. Therefore, 
we can not inspect the distributions of 33 kty and 46.3 kty exposures with 10 bins. 
Accordingly, we subdivided the whole observed zenith angle distributions separately 
into 5 bins for the respective distributions of 33 kty and 46.3 kty exposures. In table 
[|, their zenith angle distributions of various events are represented with 5 bins for 
— 1 < cos 6 < 1. Further, for comparing both the upward-directed and downward- 
directed events among them, we subdivided the data into 3 bins each for —1 < cos# < 
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and < cos 6 < 1. The region of —0.2 < cos# < 0.2 being near the horizontal direction 
was included for both distributions. 



Table 5. Chi-square tests for comparisons of the zenith angle distributions between 
33 kty and 46.3 kty exposures of e-like and fi-\ike events in the sub-GeV and the 
multi-GeV regions. 



zenith angle 


sub-GeV e-like 


sub-GeV /U-like 


multi-GeV e-like 


cos 9 


33 kty 


46.3 kty 


33 kty 


46.3 kty 


33 kty 


46.3 kty 


-1.0 0.6 


288 


349 


183 


291 


50 


41 


-0.6 0.2 


233 


344 


226 


273 


55 


73 


-0.2 ~ 0.2 


258 


306 


229 


304 


70 


95 


0.2 ~ 0.6 


226 


304 


265 


352 


75 


64 


0.6 ~ 1.0 


226 


330 


255 


410 


40 


60 


X 2 (whole) / v 
X 2 (upward) jv 
X 2 (downward) / v 


7.97/5 s.p.=0.17 
7.52/3 s.p.=0.068 
4.29/3 s.p.=0.24 


8.34/5 s.p.=0.15 
4.90/3 s.p.=0.18 
3.84/3 s.p.=0.28 


13.92/5 s.p=0.017 
5.78/3 s.p=0.13 
8.18/3 s.p=0.044 


zenith angle 


multi-GeV FC 


multi-GeV PC 


multi-GeV (FC+PC) 


-1.0 0.6 


32 


44 


32 


47 


64 


91 


-0.6 0.2 


31 


44 


44 


68 


75 


112 


-0.2 ~ 0.2 


45 


64 


91 


131 


136 


195 


0.2 ~ 0.6 


73 


78 


68 


124 


151 


202 


0.6 ~ 1.0 


49 


98 


66 


83 


115 


181 


X 2 (whole) jv 
X 2 (upward) jv 
X 2 (downward) / v 


7.11/5 s.p.=0.21 
0.015/3 s.p.=0.99 
7.10/3 s.p.=0.073 


3.93/5 s.p.=0.56 
0.33/3 s.p.=0.95 
3.64/3 s.p.=0.30 


1.36/5 s.p=0.93 
0.22/3 s.p=0.97 
1.18/3 s.p=0.76 



It is generally known that a significance probability (s.p.) of > 0.05 in a statistical 
test indicates consistency for the two data sets. In the present cases, also, s.p. values of 
little larger than 0.05 may indicate some degrees of mis-identification. In the following 
we carefully check the variances of the s.p. values of the respective distributions. In table 
[|, a small s.p. value of 0.068 was seen for the upward directed zenith angle distribution 
of sub-GeV e-like events. Further, in multi-GeV e-like events, we found an extremely 
rare s.p. value = 0.017 for the whole zenith angle distribution. The downward directed 
events were also described by a s.p.= 0.044 which was smaller than 0.05, and the upward 
directed one of s.p.= 0.13 was smaller than the other upward directed ones except for 
the sub-GeV e-like events. Although multi-GeV //-like (FC + PC) events showed very 
large s.p. values for all the distributions, the downward-directed FC events could only 
be characterized by a s.p. value of 0.073. However, all distributions of the PC events 
show large s.p. values. These varying results might be connected to the identification 
problem of muons and electrons. 

For testing our s.p. values, we used the upward through-going muon flux and the 
upward-going stopping events in — 1 < cos# < expressed with 10 bins for the former 
and 5 bins for the latter of the 79.3 kty exposure. Both these events are sure to be muons 
as they are irrelevant to the identification of muons and electrons. Similar x 2 tests were 
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done for both distributions expressed by 5 bins each from the subdivision of the 79.3 
kty exposure into the respective 33.0 kty and 46.3 kty exposures. From these data, we 
obtained s.p. values of 0.66 for 5 bins for the upward through-going muon flux, and a 
s.p. value of 0.37 for 5 bins in the upward-going stopping events. These s.p. values are 
not in conflict with being muons. From this we infer that the poor s.p. values shown in 
table |5| 5, for some of the zenith angle distributions of e-like and of /x-like events in the 
sub-GeV and multi-GeV ranges and also for the downward directed /i-like (FC) events, 
may be connected to the poor identification capability for muons and electrons in the 
33 kty and 46.3 kty exposures. Thus, we come to the conclusion that these events may 
show an appreciate mis-identification between muons and electrons. 

Accordingly, the comparison of the observed zenith angle distribution of various 
events with the expected one based on the neutrino oscillation hypothesis is still an 
open question. In addition, the comparisons of the upward-going muon flux with the 
null and the — v T neutrino oscillation hypothesis have to face the problem about 
the absolute value of the neutrino and muon fluxes at ground level for different zenith 
angles. Namely, very recently, the measurement of atmospheric vertical muon fluxes at 
ground level in the momentum ranges from 200 MeV/c to 120 GeV/c by the CAPRICE 
group Jl9| and from 0.6 to 20 GeV/c by the BESS group p0) were done by using low mass 
superconducting magnet spectrometers. The former are about 10 % to 15 % lower and 
the latter about 20 % lower than the previous experimental results. The OKAYAMA 
group |21[ using a solid iron magnet telescope has obtained the vertical muon fluxes from 
1.5 to 100 GeV/c which are consistent with both the above measurements. Also, the 
NMSU-WIZARD/CAPRICE98 |22[ balloon-borne magnetic spectrometer reports new 



measurements at several atmospheric depths in the momentum range 0.3 - 20 GeV/c. 



5. Conclusion 



Based on our MC simulation which includes fluctuations of the produced number of 
photoelectrons due to the stochastic nature of energy-loss mechanisms, we conclude 
that the total mis-identification rate for muon and electron events is larger than or 
equal to 20 % in the sub-GeV region and also at least several % in the multi-GeV 
region. In both cases the correct treatment of the particle identification would lead to 
an increase of the muon neutrino flux, namely of 20 % in the sub-GeV range and 5 % 
in the multi-GeV region. Among others, this mis-identification problem is related to 
the fact that electrons will always initiate electromagnetic cascades leading to a fuzzy 
image of the Cerenkov ring pattern. On the other hand muons can occasionally undergo 
electromagnetic interactions (knock-on production, bremsstrahlung and direct electron 
pair production). In these cases the muon track overlaps with the shower initiated by 
the secondaries created by the muon. The emerging pattern can easily be mistaken for 
an electron. 

In addition, our x 2 tests showed that the zenith angle distributions observed by 
Super K. for e-like and /i-like (FC), except /z-like (PC) events in the first half of the 
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total exposure of 79.3 kty (33.0 kty) and the second half (46.3 kty) are inconsistent 
both for the sub-GeV and multi-GeV ranges. In special, very different s.p. values 
between 0.99 and 0.073 for both observed zenith angle distributions of upward-directed 
and downward-directed /z-like (FC) events suggest that the downward going events 
suffer from mis-identifications between muons and electrons. Given these results, we 
conclude that due to the uncertainty of the identification capability for muon and 
electron events in the Super K. detector one cannot yet firmly establish the existence 
of an atmospheric neutrino anomaly. A better knowledge and precision treatment of 
the particle identification technique including discussions of the absolute intensities for 
the zenith angle dependent upward going events and the accurate atmospheric neutrino 
fluxes in general will certainly help to resolve the remaining problems. 
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